Metallic stripe in two dimensions: stability and spin-charge separation 
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The problems of charge stripe formation, spin-charge separation, and stability of the antiphase 
domain wall (ADW) associated with a stripe are addressed using an analytical approach to the t-J z 
model. We show that a metallic stripe together with its ADW is the ground state of the problem 
in the low doping regime. The stripe is described as a system of spinons and magnetically confined 
holons strongly coupled to the two dimensional (2D) spin environment with holon-spin-polaron 
elementary excitations filling a one-dimensional band. 
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In a light of experimental evidence for incommensu- 
rate spin and charge order in high-T c superconductors |l]] 
the concept of stripes, i.e. one-dimensional spatial struc- 
tures of holes condensed at an antiphase domain wall 
(ADW) of antiferromagnetic (AF) spins, has been dis- 
cussed extensively in the literature In the physics 
of high-T c superconductors magnetism plays an impor- 
tant role j[o| and it is believed that stripe formation is 
due to the strong interplay between the charge and spin 
degrees of freedom || . Various proposals for the driving 
mechanism of such a stripe phase have been made. A 
large body of them emphasize the importance of the so 
called topological doping, i.e. the breaking of the symme- 
try of the ground state in order to lower the total energy 
of the system when the charge carriers are introduced 
. For the cuprates the symmetry breaking corresponds 
to the change of the direction of the staggered magne- 
tization across the stripe, i.e. to the formation of an 
ADW. In this context the importance of the kinetic en- 
ergy of holes, both along the domain wall and transverse 
to it, has been noted in the studies of the Hubbard and 
the t-J models ^ |j. The role played by the long-range 
Coulomb field in the ordering of the stripes has been also 
discussed P|§. Numerical studies of the t-J model indi- 
cate that the (dynamic) stripe phase may be the ground 
state of the system but there is still no consensus on this 
issue pyrH]. However, the actual picture of such topolog- 
ical doping is far from being complete mainly due to the 
lack of microscopic insight into the problem. There is no 
clear understanding of the stabilization mechanism of the 
stripe as well as of the nature of low-energy excitations of 
the system. Thus it is important to study a model of AF 
spins and mobile holes within an approximation which is 
not subject to the limitations of mean field theories or 
numerical studies. 

In this paper we study analytically the problem of the 
stability of an ADW in the anisotropic t-J z model. We 
calculate the Green's function of a charge excitation at 
the ADW in a well controlled procedure. We show that 
although an empty ADW (no holes) is a highly excited 
state it is stabilized even at low doping concentrations 



due to the gain in hole kinetic energy. As a consequence 
the doped holes become spatially confined in the direc- 
tion transverse to the ADW and exhibit spin-charge sep- 
aration features. Namely, the holes behave like holons 
(charge- 1, spin-0 fermions) at the ADW and they are 
spin polarons (charge-1, spin-1/2 fermions) outside of the 
ADW. The confining potential in the transverse direc- 
tion leads to quantized states whose energy grows roughly 
~ jV 3 ■ We argue that due to the ID nature of an ADW 
holon even Trugman loops |L2j are energetically costly. 
We also compare some of our results to recent DMRG 
calculations for the isotropic t-J model || and find very 
good agreement. 

Our starting point is the t-J z model which is given by: 



n = -t 



( 5 ^ c 



Ec.) + Jj][^ 



-Ni v 
4 ' 



J J 



(1) 



where t is the kinetic energy, J is the AF exchange, 
and Ni = n,-|- + n^. All operators are defined in the 
space without double-occupancy of the sites. The natu- 
ral source of exchange anisotropy in real systems is the 
spin-orbit coupling |l3[ . Although the t-J z model is the 
strongly anisotropic limit of the more realistic SU(2) t-J 
model, it is known from numerical and analytical stud- 
ies that the low-energy physics of holes is very similar in 
both models Km . Therefore, as long as the hole dynamics 
is concerned, it should be instructive to study the t-J z 
model which captures general properties of the doped 
AFs and allows an analytical treatment of the problem. 

According to the idea of topological doping a single 
charge carrier must benefit energetically from breaking 
the symmetry of the ground state. Thus, it is natural 
to consider the problem of one hole at the ADW, in or- 
der to obtain the low-energy variables which then can be 
used for description of many-hole system. The single-hole 
problem in a homogeneous AF background is very well 
studied with analytical results and numerical data being 
in a very good agreement fll5| . The charge quasiparti- 
cle is understood as a spin polaron, i.e. a hole dressed 
by strings of spin excitations [14 15 1. In other words, 



the hole movement in a regular AF is frustrated because 
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of the tail of misaligned spins following the hole. The 
idea that an ADW can be more favorable configuration 
for holes relies on the fact that such a frustration of the 
hole's kinetic energy can be avoided for movement in- 
side the wall, such that the hole is essentially free in the 
ID structure. Then, holes populating this ID band can 
compensate the cost of magnetic energy (~ J ■ length) 
of the wall. However, this ID energy alone (~ — 2t) can- 
not overcome the energy of the spin polaron in the bulk 
E p » — 2v3£ (at t >■ J), especially if the energy cost 
for creating an ADW is taken into account. Therefore, 
the "transverse" hole movement away from the domain 
boundary must be equally important. Moreover, it is ev- 
ident that such a movement is very similar to the hole 
movement in a homogeneous AF, i.e. it also generates 
strings of misplaced spins. It is easy to see that such 
strings are also weaker, i.e. cost less energy near the 
ADW than in the bulk. Then, the coupling of the longitu- 
dinal and transverse hole movements for this many-body 
problem gives a true low-energy elementary excitation, 
which unifies the features of the ID charge carrier with 
the properties of the spin polaron. 

Let us consider a bond-centered domain wall with one 
hole injected in it. If the hole is allowed to move along the 
wall it will create a spin defect (a spinon - charge-0, spin- 
1/2 fermion, which is massive in our case because of the 
anisotropy) and then propagate freely as a holon. The 
structure shown in Fig.^ja is the natural starting point 
for the consideration of a hole whose movement along the 
stripe is free from the beginning. Fig.|l|b shows the same 
configuration schematically. It is clear that the domain 
wall not only corresponds to an antiphase shift of the 
staggered magnetization across the stripe, but also en- 
sures that the hole moves as a holon in the longitudinal 
direction. Fig^Jc shows an example of a string generated 
by the transverse hole movement. Notably, the first ele- 
ment of such a string is a spinon. These strings lead to 
an effective confinement of the hole in the direction per- 
pendicular to the stripe and to the formation of a spin- 
polaron-like cloud of spin excitations around the holon. 
This can be viewed as a generalization of ID physics to 
higher dimensions in the presence of an ADW [W . 

The bare Green's function for longitudinal hole move- 
ment in Fig. is given by 

G° XQ (k y ,Lj) = [oj-2tcos{k y ) + tO}- 1 , (2) 

so the holon band minimum is located at k y — it, where 
the index xq corresponds to the x-coordinate of the 
stripe. Then the renormalization of the Green's func- 
tion, as in the case of the spin polaron, is coming from 
the retraceable path movements of the hole away from 
ADW and back. The retraceable path approximation is 
equivalent to the self-consistent Born approximation for 
the self energy. Corrections due to terms beyond this ap- 
proximation are very small |lq ] and will be omitted. The 
full Green's function is then given by 



G Xo (k y ,oj) = [oj-2tcos{k v )-j: xo {Lu)+iO}- 1 , (3) 



where E Xo (o;) takes the form of a continued fraction 
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where is the energy of the i-th segment of the string, 
which is equal to the number of broken AF bonds (J/2 
each) associated with the segment. In the retraceable 
path approximation S Xo (w) has no k y dependence. Since 
the energy spectrum of the elementary excitations is 
given by the poles of the Green's function Eq. (||) with 
the self energy Eq. (|J), one needs to calculate and 
seek solutions of E(k y ) — 2tcos(k y ) — Yl(E(k y )) = 0. A 
standard simplification is to assume that the energy of 
the string is independent of the path of a hole and is sim- 
ply proportional to the length of the path. For the spin 
polaron this is plausible because only very few strings do 
not follow this rule. Then an analytical solution for the 
self energy is given by the ratio of Bessel functions |l5) . 
If we assume that u)\ = J/2, cjj>i = J (two broken bonds 
per segment of string, see Fig. |l|c), Eq. (|J) transforms 
to 
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with T(w) = J-u/j(r)/J- u /j-i(r), J v {r) the Bessel 
function, and r = 2^/it/J. The energy of the lowest 
pole of G{k y ,uj) Eq. (g) at J/t = 0.4 with from 
(0) versus k y is plotted in Figj| (dotted line). The zero 
energy level is chosen to be equal to the energy of a static 
hole in the configuration in Fig.|l|a. This reference to the 
t = limit is natural to show which of the magnetic 
configurations is best for optimizing the kinetic energy. 

One can consider the problem more rigorously taking 
into account the energy of each string exactly up to a 
certain length l c and applying the path-independent as- 
sumption only for I > l c . The lowest-pole energy versus 
k y for such calculations with l c — 4 is shown in Fig.^ 
(solid line) . Because of the ADW, there is less energy re- 
quired to create a spin flip near it and therefore there is 
a subset of strings having lower energy than the string of 
the same length in the bulk. It is worth mentioning that 
the energy of this ADW elementary excitation relative to 
the energy of the magnetic background is lower than the 
energy of the spin polaron in the bulk at all k y (Figj^). 
For J/t = 0.4 the gain of the energy of the holon at the 
bottom of the band over the spin polaron is about 1.5J, 
that is the energy of three broken AF bonds. Another 
informative quantity, the residue of the Green's function 
Z(k y ), is shown in Fig.|| (inset). It gives a measure of the 
amount of "bare" holon in the wave function of the ele- 
mentary excitations. One can see that a significant part 
of the initial holon at k y = n resides inside the wall and 
almost all its weight is transfered to strings at k y < ir/2. 
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Because the band is very flat at the same k y < tt/2, 
the velocity of the elementary excitations is much slower 
than the bare Fermi velocity v F — 2t sin(kp)/h. It is 
easy to show that the velocity at the Fermi level for our 
ID band is vf = v F Z(kp). For realistic values of J, t, 
and Uf around it/2 one finds Tivp ~ 50 — 100 meVA, 
which is close to the velocity of the slow mode (~ 35 
meVA) deduced from experimental data in Ref. |l7j for 
the cuprates. 

The next step is to fill the obtained ID band up to 
some Fermi momentum kp and study the stability of 
the whole system as a function of the ID hole concen- 
tration nil. Obviously, the rigid band filling neglects all 
effects of interaction between the carriers except Fermi 
repulsion. Such interactions would include attractive as 
well as repulsive terms coming from the nearest-neighbor 
hole-hole interaction, crossing of the strings, spin-flip 
exchange, and kink-kink or kink-antikink scattering of 
holons. However, as a first step towards the microscopic 
study of the metallic stripe, the rigid-band approxima- 
tions should help to establish the energy scale underpin- 
ning the stripe ground state. The total energy of the 
stripe must include the magnetic energy paid for the do- 
main wall E w — (L y — Nh)J/2 and kinetic energy of holes 
Ekin = Ylk<k F E(k y ). Here, L y is the number of sites in 
the y-direction, Nh is the number of holes, kp = imu is 
the Fermi momentum relative to 7r (band minimum), and 
nil = Nh/Ly. Thus, the total energy per hole is 

E tot /N h =_-!+- / E(k y )dk y . (6) 

Results for E to t/t versus nil at J/t = 0.4 are shown 
m Fig.| The stripe wins over the homogeneous phase 
when nil is as low as 1/4 and then it has significantly 
lower energy. Another feature of the data is that from 
nil rs 0.5 to nn = 1 the energy curve is almost flat, al- 
though it does not show a minimum. A minimum at 
nil =/= 1 will appear if repulsive short-range interactions 
are taken into account. The energy gain we obtained for 
the representative value of J/t is about 0.25i ~ 0.6 J per 
hole. The picture shown in Fig.|^ varies slowly as a func- 
tion of J/t, with the boundary of stripe stability moving 
towards smaller nn at smaller J/t. Fig|| (inset) shows 
the spatial density of holes N(x) and modulus of the 
staggered magnetization |M(x)| = |(S z (x))| in the direc- 
tion perpendicular to the stripe at the linear hole den- 
sity ny = 2/3. To extract these we calculated an explicit 
distribution of the hole and spin flips within the wave 
function of the single elementary excitation and then in- 
tegrated these quantities over k y . Our results for N(x) 
are in remarkable agreement with the same quantity for 
the single stripe obtained numerically in Ref. M for the 
case of the isotropic t - J model. 

Another issue to address is the stability of an ADW to 
the untwisting of the order parameter. In other words, 



it is necessary to show that the 7r-shift of the order 
parameter is stable against small distortions. Energy 
change of the transverse hole fluctuations by the un- 
twist is AE tr ~ +A(n\\)ip 2 , ip being the angle of the 
untwist. On the other hand, the energy of magnetic frus- 
tration caused by an ADW is decreased by the untwist, 
AE W — —B(n\\)ip 2 with the actual numbers showing that 
A > B only at large hole densities nn ~ 1. At the same 
time, the impact of the small distortion on the longitu- 
dinal hole movement is more dramatic. When the phase 
shift across the wall is different from n the hole motion 
along the wall will lead to the tail of weak spin defects 
with energy </</2 2 /4 <c J per unit length and the hole will 
be localized on the lengthscale h oc ~ (Jip 2 /t)~ 2 / 3 ^> 1. 
The corresponding energy increase due to the untwist is 
AEi ~ +C(^ 4 / 3 , which always dominates the y> 2 -term at 
small ip and ensures the local stability of the 7r-shift at 
all hole densities. 

In the case of a hole in a pure Ising background it is 
known that the hole can escape from the confinement 
potential via some high energy processes |l2| remaining 
beyond the retraceable path approximation (Trugman 
loops). However, in the case of the hole at the ADW the 
holon in the stripe has no spin, whereas the hole in the 
bulk has both spin and charge. Therefore, to leave the 
stripe and acquire a spin the hole must create a spinon 
which costs energy ~ J. In other words, a holon can only 
virtually decay into a spin-polaron and a spinon. 

In the more realistic t - J model spins are dynamic, 
which would renormalize the energies of the magnetic 
background and the holes in our picture. We believe, 
however, that the essential physics of the system will re- 
main the same. 

In conclusion, we have presented an analytical study 
of a stripe of holes at an ADW of AF spins. Longitudinal 
as well as transverse kinetic energy of holes are explicitly 
taken into account in our approach, and their role in sta- 
bilization of the stripe as a ground state of the system is 
revealed. We have provided a description of the charge 
carriers building the stripe as a system of ID elementary 
excitations, unifying the features of holons, spinons, and 
AF spin polarons. This represents a new level of under- 
standing of the structure of the stripe phase in cuprates. 
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FIG. 1. Single hole at an ADW separating AF domains. Broken AF bonds are marked by strings, (b) Same as (a). Pluses 
and minuses represent staggered magnetization, (c) String of spin flips generated by the hole. 
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FfG. 2. Energy of elementary excitation E(k y ) v.s. k y , residue of the Green's function versus k y (inset). Dotted lines are 
for path-independent string calculations, solid lines are the results of more rigorous calculations described in the text. Solid 
straight line is the energy of a spin polaron in the bulk. All energies are relative to the energy of a static hole in a corresponding 
magnetic background. 
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FIG. 3. Total energy of the system per hole versus nn . Horizontal line is the energy of free spin polarons in the homogeneous 
AF. Dotted and solid lines are the same as in Fig.^, J/t = 0.4. Inset: spatial density of holes (diamonds) and modulus of the 
staggered magnetization (circles) across the stripe at ny = 2/3, J/t = 0.35. Lines are guide to the eye. Empty diamonds are 
numerical data from Ref. 101 
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